TRANSFORMATIONS OF MARKOV PROCESSES AND CLASSIFICATION SCHEME FOR 

SOLVABLE DRIFTLESS DIFFUSIONS 



CLAUDIO ALBANESE AND ALEXEY KUZNETSOV 

Abstract. We propose a new classification scheme for diffusion processes for which the backward Kol- 
mogorov equation is solvable in analytically closed form by reduction to hypergeometric equations of the 
Gaussian or confluent type. The construction makes use of transformations of diffusion processes to elim- 
inate the drift which combine a measure change given by Doob's h-transform and a diffeomorphism. Such 
transformations have the important property of preserving analytic solvability of the process: the transition 
probability density for the driftless process can be expressed through the transition probability density of 
original process. We also make use of tools from the theory of ordinary differential equations such as Li- 
ouville transformations, canonical forms and Bose invariants. Beside recognizing all analytically solvable 
diffusion process known in the previous literature fall into this scheme and we also discover rich new families 
of analytically solvable processes. 



1. Introduction 

An "analytically solvable" Markov process can be informally defined as follows: 

Definition 1.1. A process Xt is solvable if its transition probability distribution can be expressed as an 
integral over a quadratic expression in hypergeometric functions. 

This definition is very general as it includes all the well known examples in the literature such as 
Brownian hypergeometric Brownian motion and the Omstein-Uhlenbeck, Bessel, square-root and Jacobi 
processes. It also includes a broad family of other processes which are discovered by means of the 
classification exercise in this paper. A similar classification problem but addressing the question of 
classifying all diffusion processes for which one can express the Laplace transform of the integral in 
analytically closed form, was addressed by C. Albanese and S. Lawi in [4]. 

In section 2 we briefly review the necessary definitions, facts and theorems about the diffusion pro- 
cesses. The most important objects which play a role in our constructions are Markov generators and 
transition probability densities, the speed measure and the scale and Green's functions. We also recall 
Feller's classification of boundary conditions for a diffusion process. 

In section 3 we introduce the concept of stochastic transformations as a composition of a Doob's 
h-transform and a diffeomorphism, and show how to construct a complete family of stochastic transfor- 
mations for a given diffusion process. We then show how these results can be generahzed to arbitrary 
Markov process. 

In section 4 we prove some useful properties of stochastic transformations and discuss the equivalence 
relation these transformations induce on the set of all driftless diffusions. Examples include Brownian 
and hypergeometric Brownian motions. 

In the last section 5 we generalize the concept of stochastic transformations into that of a general 
transformation of the Markov generator regarded as a second order differential operator. We introduce 
and define "Bose invariants" which are invariant under stochastic transformations and "Liouville trans- 
formations" which act on second order differential operators while preserving the Bose invariants and 
thus to obtain new families of solvable processes. We conclude by stating and proving two classification 
theorems. 

Appendix A gives necessary facts and formulas about hypergeometric functions while appendix B 
provides several useful facts about Ornstein-Uhlenbeck, square-root and Jacobi diffusions. 
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2. Background from the theory of diffusion processes 

This section is a brief introduction to the classical theory of one-dimensional diffusions: construction 
and probabilistic descriptions of analytical tools as the speed measure, the scale and Green function and 
the description of the boundary behavior of the diffusion process Xf. References on this subject include 
[10], [9], [12] and [7]. 

Let D be the (possibly infinite) interval [£'\-D^] C M, with oo < < < oo. Let Xf be a 
stationary Markov process taking values in D with transition probability function P{t, x, A) = Px{Xt € 
A). 

Definition 2.1. The probability semigroup is defined as the one-parameter family of operators 

(2.1) P{t)f{x) = I f{y)P{t, X, dy) = Eo,J{Xt) 
and the resolvent operator is defined as the Laplace transform of P{t) 

poo 

(2.2) R{X)f{x)= / e-^'P{t)f{x)dt 

Jo 

on the domain L°°{D) of bounded measurable functions / : I? — > ]R. 

We call a process conservative if P{t, x,D) = 1 for all t and all x £ D. If the process is not 
conservative, then we can enlarge the state space by adding a cemetery point A^: 

P{t,x,A^) = l-P{t,x,D). 

With this addition, the process Xt on the space D U Aqq is conservative. If / is a function on D, we will 
extend it to Z) U Aoc by letting /(Acc) = 0. 

Definition 2.2. The infinitesimal generator C of the process Xt is defined as follows: 

(2.3) Cf := ^P{0+)f = lim ^^^^^ ' ^ 

for all continuous, bounded / : D — > M, such that the limit exists in the norm. The set of all these 
functions / is the domain of C and is denoted 

Below we assume that Xt is a regular diffusion process, specified by its Markov generator 

(2.4) Cf = ]^a\x)f"{x) + b{x)r{x) - c{x)f{x) 

where the functions b{x),c{x) and a{x) are smooth and c{x) > 0, a{x) > in the interior of D. 

Every diffusion process has three basic characteristics: its speed measure m{dx), its scale function 
s{x) and its killing measure k{dx). For the diffusion specified by the generator (2.4) these characteristics 
are defined as follows: 

Definition 2.3. Speed measure and killing measure are absolutely continuous with respect to the Lebesgue 
measure (in the interior of domain D) 

m{dx) = m{x)dx, k{dx) = k{x)dx, 

and the functions m{x), k{x) and s{x) are defined as follows: 

(2.5) m(a;) = 2(7-2 (x)e-^(^), s\x) = e-^^''\ k{x) = c{x)m{x) = 2c{x)(j-^{x)e^^''^ 
where B(a;) := 2a~'^{y)b{y)dy. 

Remark 2.4. We denote by m{x) a density "^l^^ ■ The same applies to the kUUng measure k{dx). 
The functions m, s and k have the following probabilistic interpretations: 
• Assume fc = 0. Let := inf{t : Xt = z} and (a, b) C D. Then 

s(b) — s(a) 

We say that Xt is in natural scale if s{x) = x. In this case (if the process is conservative) Xt is 
a local martingale. 
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• The speed measure is characterized by the property according to which for every t > and 
a; e -D, the transition function P{t, x, dy) is absolutely continuous with respect to m{dy): 

P{t,x,A) = j p{t,x,y)m{dy) 

A 

and the density p{t, x,y)is positive, jointly continuous in all variables and symmetric: p{t, x,y) = 
p{t, y, x). Notice that the transition probabiHty density p{t, x, y) is the kernel of the operator 
P{t) with respect to the measure m{dy). 

• The killing measure is associated to the distribution of the location of the process at its Ufetime 
C := ini{t : Xt ^ D): 

t 



V^{Xc- €A\C<t) = jdsj p{s,x,y)k{dy). 



From this point onwards we assume that there is no killing in the interior of domain D, or namely that 
c{x) = 0. 

Remark 2.5. Notice that the scale function can be characterized as a solution to equation 

£s{x) = 0, 

and s'{x) is proportional to the Wronskian W^^^^^ {x), where </9i, <^2 are any two linearly independent 
solutions to 

Let r be the stopping time with respect to the filtration {J^t}- The process X^/^f is called the process 
stopped at r and is denoted by X^. 

The following lemma is required in the next sections: 

Lemma 2.6. Let T = inf{t > : Xt ^ int(£')} be the first time the process Xt hits the boundary of 
D. Then for each x G D, = s{Xj) is a continuous V ^-local martingale. 

We refer to [12], vol. 11, p.276 for the proof of this lemma . 

The speed measure and the scale function are defined in terms of the coefficients of the generator C. 
Also the inverse of the above claim holds true: if the generator C of the process Xt can be expressed as 

(2.6) Cf=^^ff^=DM, 

then the speed measure and the scale function define the generator of the process Xt (and thus determine 
the behavior of Xt up to the first time it hits the boundary of the interval). The boundary behavior of the 
process Xt is described by the following classical result (see [9J,[10]): 

Lemma 2.7. Feller classification of boundary points. Let d G {D^,D^). Define functions 

R{x) = m{{d, x))s'{x) and Q{x) = s{x)m{x). Fix small e > (such that + e G D). Then the 

endpoint is said to be: 

'regular if Q e L^{D\D^ + e), R e {D^ , + e) 

exit if Q i L^{D^,D^ + e), R e L^D\ + e) 

entrance if Q G L^{D^ , + e), R^ L\D\D^ +e) 

mturalif Q ^ L\D\D^ + e), R^L\D\D^ +e) 

The same holds true for D"^. 

Next we elaborate on the probabilistic meaning of different types of boundaries. 

Regular or exit boundaries are called accessible, while entrance and natural boundaries are called 

inaccessible. 

An exit boundary can be reached from any interior point of D with positive probability. However it is 
not possible to start the process from an exit boundary. 



(2.7) 
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The process cannot reach an entrance boundary from any interior point of D, but it is possible to start 
the process at an entrance boundary. 

A natural boundary cannot be reached in finite time and it is impossible to start a process from the 
natural boundary. The natural boundary Di is called attractive if Xi as t ^ oo. 

A regular boundary is also called non-singular. A diffusion reaches a non-singular boundary with 
positive probabiUty. In this case the characteristics of the process do not determine the process uniquely 
and one has to specify boundary conditions at each non-singular boimdary point: if m({D*}) < oo, 
fc({Z)*}) < oo, then the boundary conditions are 



'g{D^)m{{D^}) - ^ + f{D^)k{{D^}) = 0, 
a{D^)m{{D'}) + ^ + f{D^)k{{D^}) = 0. 



(2.8) 

ds{x) 

where g := £/ for / e V{C). 

The following terminology is used: the left endpoint is called 

• reflecting, if m{{D^}) = k{{D^}) = 0, 

• sticky, if to({L>1}) > 0, k{{D^}) = 0, 

• elastic, if to({£)1}) = 0, ki{D^}) > 0. 

A diffusion process X spends no time and does not die at a reflecting boundary point. X does not 
die, but spends a positive amount of time at a sticky point (which in the case m{{D^}) = oo is called 
an absorbing boundary - the process stays at forever after hitting it). X does not spend any time 
at elastic boundary - it is either reflected or dies with positive probability after hitting (in the limit 
k{{D^}) = oowe call a killing boundary, since that X is killed immediately if it hits D^). 

Let the interval D be an infinite interval, for example of the form [D^,oo). We say that the process 
Xi explodes if the boundary = oo is an accessible boundary. Using the previous lemma one can see 
that the process explodes if and only if for some e > 

(2.9) R{x) = m{{D^ + e,x))s'{x) e L\D^ + e,oo). 

In section 3 below, we construct two linearly independent solutions to the ODE 

(2.10) jCip{x) = X^fix), A > 0, xGD. 

The probabilistic description of these solutions is given by the following lennma (see [12], vol. II, p. 
292): 

Lemma 2.8. For A > there exist an increasing f'i^ (x) and a decreasing ip~^ (x) solutions to equation 
(2.10). These solutions are convex, finite in the interior of the domain D and are related to the Laplace 
transform of the first hitting time as follows: 



(2.11) E,{e-^"') = 



x < z, 
x > z. 



The functions <^J(a;) and ip^ (x) are also called the fundamental solutions of equation (2.10). These 
functions are linearly independent and their Wronskian can be computed as follows: 

(2.12) W^.t.-(-) = '^^A (-) - (-)'^ = -A.'(x), 

thus the Wronskian with respect to Dg = d/ds{x) is constant: 

The following theorem due to W. Feller characterizes boundaries in terms of solutions to the equation 
(2.10) and will be used in section 3: 

Theorem 2.9. (i) The boundary point is regular if and only if there exist two positive, decreas- 
ing solutions (fi and <p2 of (2.10) satisfying 

(2.14) limv5i(x) = 0, lim ^^i^ = -l, lim ^2{x) = I, lim = 0. 

x^D^ ^ x^D^ ds{x) x^D^ ^ ' x^D^ ds{x) 
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(ii) The boundary point is exit if and only if every solution of (2.10) is bounded and every positive 
decreasing solution satisfies 

(2.15) lim ^i(a;)=0, lim < 0. 

x^D'^ x^D^ ds[X) 

(iii) The boundary point is entrance if and only if there exists a positive decreasing solution (fi 
of (2.10) satisfying 

(2.16) lim <pi(x) = l, lim '^pM = q, 

x^D^ x->£)2 ds{x) 

and every solution of (2.10) independent of ipi is unbounded at D'^. In this case no nonzero 
solution tends to as x ^ D^. 

(iv) The boundary point is natural if and only if there exists a positive decreasing solution ipi of 
(2.10) satisfying 

(2.17) lim ipAx)=0, lim = 0, 

x^D^ x->£)2 ds{x) 

and every solution of (2.10) independent ofipi is unbounded at D"^. 
In cases (i) and (ii), all solutions of (2.10) are bounded near D2 and there is a positive increasing 
solution z such that lim z{x) = 1. In cases (iii) and(iv) every positive, increasing solution z satisfies 

lim zix) = 00. 

x^D^ 

Another important characteristic of a Markov process is the Green function: 

Deiuiition 2.10. The Green function G{X, x, y) is defined as the Laplace transform of p{t, x, y) in time 
variable: 

00 

(2.18) G{\x,y):= J e-^'p(t,x,y)dt. 



The Green function is symmetric and it is the kernel of the resolvent operator R{X) = (£ — A)~^ with 
respect to m{dx). 



The Green function can be conveniently expressed in terms of functions (x) and ip^^ (x) as: 
(2.19) G{X,x,y) = 



wx^vt{y)vl{x), y<x. 



A diffusion is said to be recurrent if Vx{Hy < 00) = 1 for all x,y G D. A diffusion which is 
not recurrent is called transient. A recurrent diffusion is called null recurrent if Ex{Hy) = 00 for all 
x,y G D and positively recurrent if Ex{Hy) < 00 for all a;, y e .D. The following is a list of useful facts 
concerning recurrence, Green function and speed measure: 

• Xt is recurrent if and only if lim G{\, x, y) = 00. 

• Xt is transient if and only if lim G(A, x, y) < 00. 

• is positively recurrent if and only if /^(Z)) < 00. In the recurrent case we have: lim XG{X,x,y) = 
^^j^ and the speed measure m{dx) is a stationary (invariant) measure of X^: 

mP{t){A):= / m{dx)P{t,x,A)=m{A). 
J A 
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3. Stochastic transformations 

Let's ask a question: how can we transform a Markov diffusion process Xi in such a way that the 
transition probability density of the transformed process can be expressed through the probability density 
of Xtl First of all let's take a look at what transformations are available. 

Given a stochastic process {Xt , P) one has the following obvious choice: 

• change of variables (change of state space for the process): 

x, = (x„p)^r, = (y(xo,P), 

where Y{x) is a diffeomorphism Y : ~^ 

• change of measure: 

(Xt,P)^(X(,Q), 

where measure Q is absolutely continuous with respect to P: rfQj = ZtdPt- 

• time change: 

Xt = {Xt,J^uP) ^Xt = {Xr„J^r,,P), 

where tj is an increasing stochastic process, tq = 0. 

• and at last we can combine the above transformations in any order. 

In this section we will discuss only the first two types of transformations. The stochastic time change 
can be efficiently used to add jumps or "stochastic volatility" to the processes and is important for appli- 
cations to Mathematical Finance (see [3], [2]). 

The next obvious question is: do these transformations preserve the solvabihty of the process? The 
answer is always "yes" for the change of variables transformation (we assume that the function Y{x) is 
invertible): the probability density of the process Yt = Y{Xt) is given by: 

PY{t,yo,yi) = Px{t,X{yo),X{yi)), yi e Dy, 

where X = X{y) = Y~^{y) and the speed measure of Yt is 

mridy) = mY{y)dy = mx{X{y))X'{y)dy. 

What can we say about the measure change transformation? We want the new measure Q to be 
absolutely continuous with respect to P, thus there exists a nonnegative process Zt, such that 

The process Zt must be a (local) martingale. Under what conditions on Zt does this transformation 

preserves solvability? 

Informally speaking, the probability density of the process X^ = (Xj, Q) is given by 

pjs:Q(t,a;o,a;i)mxQ(.'ri) = E'^{5{xf_^^ - xi)\Xs = xo), 

where we used the fact that the transformed process is stationary. Using a formula of change of measure 
under conditional expectation we obtain 

Px<i{i,xo,xi)mxQ{xi) = ^E^{Zs+t5{Xfj^t - xi)\Xs = .tq). 

Since we want the transformed process to be Markov, we see that this can be the case if and only if Z^ 
depends only on the value of Xg, thus the process Zt can be represented as 

Zt = h{Xt,t), 

for some positive function h{x,t). 

The next step is to note that the dynamics of Xt imder the new measure is: 

dXt = dt + a{Xt)dW^^. 
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Now we see that if we want the transformed process to be stationary, hx/h must be independent of t, 
which gives us the following expression for the function h{x, t): 

h{x, t) = h{x)g{t). 

Since we don't want to introduce killing in the interior of the domain D, the process Zf = h{Xt)g{t) 
must be a martingale (at least a local one), thus we have the following equation 

'^h{x)+g{t)Cxh{x) = Q. 

After separating the variables we find that there exists a constant p such that g{t) = e"''*, and the function 
h{x) is a solution to the following eigenfunction equation: 

(3.1) Cxh{x) = ph{x). 
This discussion leads us to our main definition: 

Definition 3.1. The stochastic transformation is a triple 

(3.2) {p,h,Y} 

where p and h{x) define the absolutely continuous measure change through the formula 

(3.3) dP? = exp{-pt)h{Xt)dPu 

and Y{x) is a diffeomorphism Y : ^ Dy CR, such that the process (Y^, P'*) = is a 

driftless process. 

Remark 3.2. We deliberately do not require Yt to be a (local) martingale, since as we will see later Yt is 
not conservative in general. Though using lemma (2.6) we see that Yt is a driftless process if and only if 
Yf stopped at the boundary of Dy is a (local) martingale. 

3.1. Main Theorem. The following theorem gives an exphcit way to find all stochastic transformations 
defined in (3.1): 

Tlieorem 3.3. Let Xt be a stationary Markov diffusion process under the measure P on the domain 
Dx C M and admitting a Markov generator Lx of the form 

Assume p > 0. 

Then {p, h, Y } is a stochastic transformation if and only if 
ih{x) = ciip+{x) + C2(p-{x), 

where Cj e M are parameters ,ci, C2 > 0, C1C4 — C2C3 ^ 0, and ^p'^{x) and ^p{x) are increasing and 
decreasing solutions, respectively, to the differential equation 

Cxv{x) = pip{x). 
Before we give a proof of this theorem we need to present some theory. 

3.2. Doob's h-transform. 

Definition 3.4. A positive function h{x) is called p-excessive for the process X if the following two 
statements hold true: 

(i) e-P'E{h{Xt)\Xo = x)<h{x) 

(ii) \iraE{h{Xt)\Xo = x) = h{x). 

An p-excessive function h is called p-invariant if for a&x & Dx and t>Q 

e-P*E{h{Xt)\Xo = x) = h{x). 

Remark 3.5. Function h{x) is p-excessive (p-invariant) if and only if the process cxi'){—pt)h{Xt) is a 
positive supermartingale (martingale). Thus, if a function h{x) is zero at some xq, then /i = in Dx- 



(3.5) 



g 
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One can check that for every xi G functions (p+(a;), fp{x) and x Gx{p,x,xi) are p- 
excessive. As the following lemma shows, these functions are minimal in the sense that any other exces- 
sive function (except the trivial example h = C) can be expressed as a linear combination of them (see 
[7]): 

Lemma 3.6. Let h{x) be a positive function on such that h(xo) = 1. Then h is p-excessive if and 
only if there exists a probability measure v on = [D^ , D'^], such that for all x € 



^£,2] Gx{p,xq,xi) 

Gx{p,x,xi) , , fpjx) ^ v?+(x) 



Measure v is called the representing measure ofh. 

Deftnition 3.7. The coordinate process Xt under the measure P'' defined by 

(3.6) P\A\Xo = x) = E [e-P'^^^MX^ = ^ , 

is called Doob's h-transform or p-excessive transform of X. We will denote this process by {X, P'^) (or 
in short X^). 

The following lemma gives the expression of the main characteristics of the /i-transform of X: 

Lemma 3.8. The process X^ is a regular diffusion process with: 

• Generator 

(3.7) Cxh = j^^xh - p 
with drift and diffusion terms given by 

hxix) 

(3.8) hxh{x) = hx{x)+a\{x)^^, axh{x) = ax{x). 

• The speed measure and the scale function of the process X^ are 

(3.9) mx>^{x) = h^{x)mx{x), s'^^^x) = h-^{x)s'{x). 

• The transition probability density (with respect to the speed measure ruxh {dxi)) 

(3.10) px>'{t,xo,xi) = — rpx{t,xo,xi). 

h{xo)h[xi) 



The Green function 



(3.11) Gxh{\,XQ,xi) = f^(^^^^^^^^ ^x{p+ \xo,xi). 
• The killing measure 

I. (^^\ mx{x)v{dx) 

(3.12) kxh[dx) = - — — — r. 

Gxh(0,Xo,x} 

Proof. We will briefly sketch the proof. Using the formula (3.6) we find the semigroup for the process 

Pxn{t)I{x) = e-P'^^Px{t){hf){x), 

and from this formula we find the expression for the Markov generator and formulas for the drift and 
diffusion. Now, 

J (^x>^{y) J (^xiv) 
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from which formula we compute the expressions for the speed measure, scale function, probability den- 
sity and Green function. 

Finally, let's prove the formula for the killing measure: the infinitesimal killing rate Cxh can be 
computed as Cxh 1, thus we find: 

Cxh! = rCxh - p = -Cx / — Mdxi) - p. 

Now assuming that ^{dx) = v{x)dx and using the fact that the Green function satisfies the following 
equation 

jCxGx{p,x,xi) = pGx{p,x,xi) + d{xi - x), 

we find 

If pGx{p,x,xi) + S{xi-x) v{x) 

cxh{x) = Lxhl = - I -r- — ^ ^ v{xi)dxi - p= — r, 

hJ[D\D^] Gx{p, Xo,xi) h{x)Gx{p,Xo,x) 

thus the kilUng measure can be computed as 

7 ^ ^ ^ ^ ^^ mxix)v{dx) 

kxh{dx) = cxh{x)mxh{x)dx = - — — — r, 

Gxh[U,Xo,x) 

which ends the proof. □ 

Remark 3.9. Note that we have a nonzero killing measure on the boundary of Dr,. and no killing in the 
interior of if and only if the representing measure i'{dx) is supported on the boimdary of D^, that is 

h{x) = u{{D'})ip-{x) + u{{D^})^+{x) = ci<^;(a;) + c^cp+ix). 

Remark 3.10. To understand the probabilistic meaning of the Doob's h-transform it is useful to consider 
the procedure of constructing new process by conditioning Xf on some event (the event A can be that 
the process stays in some interval or that it has some particular maximum or minimum value). The 
mathematical description follows: 

Let the probability density Pt{x, y\A{to, ti)) be the conditional density defined as follows: 

Pt{x,y\A{to,h))dy = P{Xt+st e dy\Xt = x,A{to,ti)). 

Assume that the event A{tQ, ti) is JF^^ -measurable and it satisfies the semigroup property: P(^(<:o, ti)\B) ■ 
P{A{to,t') f] A{t' ,ti)\B) for any event 5 and f G (io,ii). 
Let the probability Tr{x, t; A{to, ti) be defined as 

Tr{x,t;A{to,ti)) = P{A{to,ti)\X{t) = x). 

Then function tt satisfies the following backward Kolmogorov equation: 

dn d7r{x,t;A{t,T)) 1 ^ d^7r{x,t; A{t,T)) dTT{x,t; A{t,T)) 

9t + ^^'^ = dt + 2^^ dx-^ + ^(^^ dx = 

The boundary conditions depend on event A{t, T). 

The conditioned drift b{x; A{t, T) and conditioned volatility (t{x; A{t, T) are given by 

2, Tix{x,t;A{t,T)) 



(3.13) h{x; A{t,T)) = b{x) + (j\x) 



7T{x,t;A{t,T)) ' 
a{x;A{t,T)) = a{x). 

Lets consider some examples of the h-transform: 

(i) Brownian bridge through /i-transform Let Xf = Wt be the brownian motion. Fix some xq and 
consider the event A{t; T) = {Xt = xi}. Then the function 7t{x, t; A{t, T)) is the probability 
density of Brownian motion: 

1 / (.Tl - xY 

Tr{x, t; A{t, T)) = pT-t{x, xi) = , ^ exp 
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One can show that X'^ is just the Brownian bridge - brownian motion conditioned on the event 
Xt = x\ (note that this process is not time-homogeneous since in this case we are not using the 
p-excessive transform). 

The conditional drift of the process X^, computed by the formula (3.13) is equal to 

TT 1 — t 

which is another way to prove that X^ is a Brownian bridge. 

(ii) Brownian Motion, Xt = Wt. Increasing and decreasing solutions to equation 

are given by 

Note that in this case the process e~''*h{Xt) is a martingale. Let h{x) = iPp{x), then the process 
Xj^ = Wf^ + \/2pt is Brownian motion with drift. Note that since e~P*h{Xt) is a martingale, the 
transformed process is still conservative, but it has a completely different behavior: for example 
the X^ ^ +00 as f — > oo and D"^ = 00 becomes an attractive boundary. This is a typical 
situation for the p-excessive transforms: as we will see later, the following result is true: the 
transformed process X'^ is either nonconservative with killing at the boundary, or in the case it 
is conservative, Xf converges to one of the two boundary points as t — > 00. 
We will return to the example of Brownian motion later in section (4. 1). 

(iii) Let Xt be a transient process with zero kilhng measure. Then Xt or Xt with 
probability 1 as f — > C- Let X+ be the iPo{x) transform of X. Note that iPo{x) is a constant 
multiple of Px{Xt D"^ as t ^ C), thus X+ is identical in law to X, given that Xt — > D"^ as 
t ^ C,,OT Otherwise X+ has the property: 

P,(limX+ = D2) = i. 
t 

3.3. Proof of the main tlieorem. Now we are ready to give the proof of the main theorem (3.3): 
Proof. Lets prove first that the process 

Zt = e-P'h{Xt) = e-P'id^+ix) + C2'fi-{x)) 
is a positive supermartingale. By applying Ito formula we find: 

dZt = {~pZt + e-P\Cxh){Xt)) dt + e-P'h'{Xt)a{Xt)dWt = erP^' {X,)a{Xt)dWt, 

since Cxh = ph. Thus Zt is a local martingale. Since it is also positive it is actually a supermartingale 
(by Fatou lemma). Thus 

(3.14) E{e-P%{Xt)\Xo = x) < h{x), 

and ex.p{—pt)h{Xt) correctly defines an absolutely continuous measure change. 

Lemmas (3.6) and (3.8) show that the converse statement is also true: if function h{x) can be used 
to define an absolutely continuous measure change, then eP*h{Xt) is a supermartingale, thus h{x) is a 
p-excessive function. Since we want the transformed process X'^ to have no killing in the interior of D^, 
representing measure ^{dx) must be supported at the boundaries (see remark (3.9)), thus we have the 
representation h{x) = ciipp{x) + C2iPp{x). 

To prove the second statement of the theorem, we note first that Y{Xt) has a zero drift if and only if 
it is in the natural scale: thus we need to check that function Y{x) given by formula (3.5) is equal to the 
scale function Sx^ (x) (up to an affine transformation), which we can check by direct computation 

_ d^f _ 'p'{x)h{x) - h'{x)tp{x) _ _ 

^ ' " dx [h{x) J ~ h^{x) ~ h^x) ~ ^' 

where (p is an arbitrary solution of Cx'fi = PV linearly independent of h, (thus it can be represented as 
ifi = c^ip'^ + C4^~ with C1C4 — C2C3 ^ 0). This ends the proof of the main theorem. □ 
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Remark 3.11. The statement that Y{Xi;) is a driftless process is an analog of the following lemma: 

Lemma 3.12. Let Q = P, and Zt — ^p^. An adapted cadlag process Mt is a P-local martingale if and 

only if Mt/ Zt is a Q-local martingale. 

Note that this lemma does not assume that the process is a diffusion process and one could use it to 
define stochastic transformations for arbitrary Markov processes. One could argue as follows: e~^*h{Xt) 
and e^''*(c3!^+(Xt) + c4ip~{Xt)) are local martingales under P, thus if we define a measure change 
density Zt = e~P*^h{Xt) the process 



Y{Xt) = 



e-P^c^^+jXt) + c^^-jXt)) _ C3^+{Xt) + c,^-{Xt) 



Zt h{Xt) 

is a Q local martingale. However, as we will see later, in some cases Zg is not a martingale, thus the 
measure change is not equivalent in general and we can't use this argument. 

3 .4. Generalization of stochastic transformations to arbitrary multidimensional Markov processes. 

The method of constructing stochastic transformations described in theorem 3.3 can be generalized to 
jump processes and multidimensional Markov processes: let Xt be a stationary process on the domain 
Dx C M"^ with Markov generator Cx- The first step is to find two hnearly independent solutions to the 
equation 

CxV^ pip- 
Assume that for some choice of Ci, C2 the function h{x) = Ci(fii{x) + C2'P2{x) is positive. Then one has 
to prove that the process 

Zt = e-P'h{Xt) 

is a local martingale (or a supermartingale), thus it can be used to define a new measure P'* by the formula 
(3.6). Then we define the function Y{x) : K'' ^ R as 

h{x) 

Now one can check that the generator of Xj^ is given by 

^xi^ = \^xh - p, 

and to prove that the (one-dimensional) process Y{Xt) is driftless one could argue as follows: 

Cx.Yix) = [ICxh -p)l = \cx {hi) - pI ^ \cx,> - = 0, 

since (p{x) = c-i,ipi{x) + C4 1^2(2;) is also a solution to Cx<p = p<p>- 

Note that this "proof" does not use any information about the process Xt and thus it is very general. 
For example in [2], [5], [6] authors use this method to construct solvable driftless processes on the lattice. 

4. Properties of stochastic transformations and examples 
The following lemma summarizes the main characteristics of the process Yt. 
Lemma 4.1. The process Yt = {Y{Xt) , P'') is a regular diffusion process with 



• Generator 



1 2, .d^ 



(4-1) ^^=o^r(2/) , 2 



where volatility function is given by: 
(4.2) <Jy{Y{x)) = ax{x)Y'{x) = ^xW^^^ 

(W{x) = s'x{x) is the Wronskian oftf^, ip~). 
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• The speed measure and the scale function of the process are 

(4.3) mriy) = = V- 

• The transition probability density (with respect to the speed measure my {dyi)) is 

(4.4) pY{t,yo,yi) = pxh{t,xo,xi) = ^^^^^^^^^^ px(^, a^o, a^i), 

where xji = Y{xi). 

• The Green function 

(4.5) Gr(A,2/o,2/i) = Gxh{\,xo,xi) = — — ^ — ^Gx{p + \xq,Xi). 

n{xo)n{xi) 

Lemma 4.2. Let Xt be a diffusion process on = [D^, D"^] with both boundaries being inaccessible. 
If 02 = (ci = 0) the domain Dy of the process Yt is an interval of the form [yo, oo) or (— oo, y^], where 
Vo = cs/ci {yo = C4/C2). In the case ci ^ and C2 7^ 0, the domain Dy is a bounded interval of the 
form [yo,yi] = [£4/02,01/03] or [yo,yi] = [01/03,04/02]. 

Proof. Remember that 



Yix) 



ci<^;J'(a;) + C2(ppix) 



The statements of the lemma follows easily from the fact that is a positive increasing function, finite 
at and infinite at D'^ and ip~ is positive decreasing function, infinite at and finite at D"^ (see 
2.9). □ 

Lemma 4.3. The processes X^ and Y are transient. 

Proof. This follows from the fact that 

lim G^.(A,a.o,xi) = ^^^-2_Gx{p,xo,x,) < 00, 

thus X^ is a transient process (see section 2). □ 

The previous lemma tells us that with probability one X^ (and Y-^) visits every point in its domain only 
a finite number of times. Thus it converges ast ^ 00. Since with probabiUty one it can not converge to a 
point in the interior of D^, it must converge to the point on the boundary (or to the cemetery point if 
the process is not conservative). Actually an even stronger result can be obtained by means of martingale 
theory: 

Lemma 4.4. If the process Y^ is conservative, then 

(4.6) p''(iim yt = yoo) = i. 

and Voo is integrable. The random variable Y^ is supported at the boundary of the interval Dy: in the 
case Dy = [yi , 00) we have 1^ = y\ a.s., while in the general case Dy = [yi , 2/2] we have that Yt 
converges to Y^ also in Li and the distribution ofY^ is 

(4.7) P(y^=y2|Fo = 2/o) = ^^^, P(F^=yi|Fo = yo)= 

y2 - 2/1 2/2 - yi 

Proof. If Yt is conservative, then If is a local martingale bounded from below (or above), thus it is a 
supermartingale bounded from below (or a submartingale bounded from above), thus it converges as 
t 00. If Dy = [2/1,2/2] and the process Yt is bounded, then it is a uniformly integrable martingale and 
it converges to Y^ also in Li, moreover 

Yt = E{Y^\J^t), 

thus 2/0 = EYt = EYoo, from which we can find the distribution of Y^o. □ 

Lemma 4.5. The stochastic transformation {X,P) (Y,P'*) given by {p,h{x),Y{x)} is invertible. 
The inverse transformation {Y, P^) t-^ {X, P) is given by {—p,l/h{x),X{y)}. 
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Proof. If the process Xt is conservative with inaccessible boundaries, the process will be a transient 
driftless process (possibly with accessible boundaries). Let's check that the process 



is a P martingale: 



since the initial process Xt is conservative and E^l = 1. Thus the transformation {—p, l/h{x),X{y)} 
maps the process Yt = {Y{Xt),P^) back into the process Xt = {Xt, P). □ 

Definition 4.6. We will say that {X, P) and {Y, Q) are related by a stochastic transformation and will 
denote it by writing 

X ^Y 

if there exists a stochastic transformation {p, h, Y} which maps {X, P) (Y, P'^). 
Lemma 4.7. The relation X ^ Y is an equivalence relation. 

Proof. We need to check that for all X, Y and Z 

(i) X ~ X 

(ii) if X ~ y, then y ~ X 

(iii) if X - y and y - Z, then X ^ Z. 

The first property is obvious; the second was proved in the previous lemma. To check the third, let 
{p\,h\{x),Y{x)} {{p2,h2{y),Z{y)}) be the stochastic transformation relating X and Y (Y and Z). 
Then {p\ + p2,h\ {x)h2{Y{x)), Z{Y{x))} is a stochastic transformation mapping X ^ Z. □ 

Thus relation divides all the Markov stationary driftless diffusions into equivalence classes, which 
will be denoted by 

(4.8) aji(x) = $m(x,p) = {(y,Q): (y,Q) ~ (x,p)}. 

Later in lenmia 5.7 we give a convenient criteria to determine whether two processes X and Y are in 
the same equivalence class (can be mapped into one another by a stochastic transformation). 

Now we illustrate with some real examples the usefulness of the concept of stochastic transformation. 
We will review some well known examples (geometric Brownian motion, quadratic volatility family, 
CEV processes) and show how these processes can be obtained by a stochastic transformation, and in 
the case of Ornstein-Uhlenbeck, CIR and Jacobi processes we will construct new famihes of solvable 
driftless processes and study their boundary behavior. 

4.L Brownian Motion. Let Xt = Wt be the Brownian motion process. Then the Markov generator is 

2dx'^' 

Fix p > 0. Functions and (p~ are given by: 

(4.9) (p+(a;) = e^% = e"^- 

The next step is to fix any two positive ci, C2 and set h = ci<^+ + C2ifip . Note that e~''^h{Xt) is a 
martingale (a sum of two geometric brownian motions). 

We consider separately two cases - (a) one of ci , C2 is zero {Dy is unbounded) and (b) both ci , C2 are 
positive. 

Let's assume first that ci = and C2 = 1, thus 

h{x) = ^-{x) = e-^P^. 

Thus the fimction Y{x) is 

y(x)=^i^^4±^l^ = C3e^v^^ + C4, 
h[x) 
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and its inverse x = X(y) is 



x{y) = TTi^ log 



2y/Tp |C3| 

Using formula 4.2 we find the volatility function of the process 1*: 

(4.10) aY{y) = <Tx{X{y))^^^^C{y-y^), 

and — 2/1 is the well known geometric brownian motion. 

Note that this process is a martingale, it is transient and lim Yt = y\ a.s. 

Now let's consider the general case: ci > 0, C2 > 0. Assume that C4/C2 > c^/ci. Then 

CsV'p (a;) + Ciip-{x) _ cge^v^^ + C4 



^ ^ ci<^p (x) + C2<Pp (a;) cie^v^x + C2 ' 



X{y) = (log f ^ ) + log (IC4/C2 - 2/1) - log (I2/ - C3/C1I) 



and 



2^ V V|ci 
Thus the derivative X'{y) is 

1 



C(c4/C2-y)(2/-C3/ci)' 

and again using formula (4.2) we find the volatility function of the process Yt 

(4.11) <yY{y) = C{ci/c2-y){y-cs/ci) = C{y2-y){y-yi), 2/2 > 2/1 • 

and we obtain the quadratic volatility family. In this case the process Yt is a uniformly integrable martin- 
gale. As t ^ 00 the process Yi converges to the random variable Fqo with distribution supported on the 
boundaries 2/1, 2/2 and given by equation (4.7). 

Remark 4.8. We have proved so far that starting from Brownian motion we can obtain the martingale 
processes with volatiUty function 

aviy) = 022/^ + ai2/ + oo, 

where the polynomial n2?y^ + aiy + 0,0 is either linear or has two real zeros. However theorem 5.7 tells 
us that the process Yt is related to Wt for all choices of coefficients Oj. Let's consider the example of the 
volatiUty function 

(4.12) c7y(2/) = 1 + 2/2 

to understand what happens in this case. 

The process Yt with volatility (4.12) is supported on the whole real line M and one can prove using 
Khasminskii's explosion test that this process explodes in finite time. We can obtain this process starting 
from Brownian motion by an analog of the stochastic transformation with p < 0. For example, let 
p = — I . Then solutions to equation 

- , , 1 cP(p(x) , , 

are given by functions sin(a;) and cos(a;). Let's choose 

h{x) = cos(a;), Y{x) = ^j^^^'^ = tan(a;). 

Then Y'{x) = ^^^^^-^ = 1 + tan^(a;) — 1 + y"^, and we obtain the volatility function given in (4.12). 

Function Y{x) = tan(a;) is infinite when x = n/2 + kir, k G Z, thus the process Yj = {Y{Wt),P^) 
explodes in finite time. 

Thus we have proved the following 

Lemma 4.9. Starting with Brownian motion Wt we can obtain the class of quadratic volatility models: 

(4.13) m{Wt) = {Yt : ayiy) = a^y^ + a^y + ao}. 
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4.2. Bessel processes. Let X^ be a Bessel process defined by generator 

dx 2 dx"^ 

We assume that a = || — 1 > (thus the process never reaches the boundary point x = 0). Then we 
prove the following: 

Lemma 4.10. As a particular case one can cover the CEV (constant-elasticity-of-variance) models with 
volatility function ayiu) = c{y — yo)^- 

Proof. Choose p = 0. Then the eigenvalue equation is 

dip(x) 1 2 (P^(x) , , 
a , ^ + -a^x ,\ ' = pipix) = 0. 
dx 2 dx^ '^^^ ' 

The two linearly independent solutions are <^o (^) — ^ Vo (^) — x~". We put C2 = 0, thus 

h{x) = a;~" and Y{x) is 

Y(x) = = A + Bx . 

Thus we can express x = X{y) = ci(y — j/o) ° and we find that 

n^(y))-(X'(?y))-i = C2(y-yo)^-i. 
Note that since a > 1, the power 1 — ^is positive. Using formula (4.2) we can compute the volatihty 

oyin) = a^/X(^Y'iX{y)) = c{y - yof , 
where 61= 1- ii. □ 

Remark 4.1 1. Note that we were able to compute the explicit form in these cases because that we could 
find the inverse function x = X{y) explicitly. This is not the case for most apphcations, though, and we 
often have to use numerical inversion instead. 

4.3. Omstein-Uhlenbeck processes - OU family of martingales. Let Xt be the Ornstein-Uhlenbeck 
process 

dXt = {a- bXt)dt + adWt, 

discussed in detail in section B.l. 

Without loss of generality we can assume that a = (otherwise we can consider the process Xt — |). 
Process X-^ satisfies the following property, Xt has the same distribution as —Xf. 

(4.14) P,{Xt e A) = P-^{-Xt e A). 

Functions ip~ {x) and (a;) are solutions to the ODE: 



One can check that the function (p^ {x) is given by 



r(i + |;) r(jy ; 

and due to the symmetry of Xt (see (4.14)) we have = (p~{—x). 

Remark 4. 12. To prove formula (4. 16) one would start with the Kummer's equation (A.l 1) for M( ^ , ^,z) 
and by the change of variables z = ^x"^ reduce this equation to the form (4. 15). 
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We see that for x > the function ip~{x) is just U{^, -^x^), where U is the second solution to 
the Kummer's differential equation (A.13). As we see in the next section, functions M and U are related 
to the CIR process , since the square of Omstein-Uhlenbeck Yt = is a particular case of CIR process: 

dYt = {a^ - 2bYt)dt + 2a^/YtdWt. 
The asymptotics of <^p (a;) as a; — > oo can be found using formula (A.15): 

(4.17) <^^(a;) ~ Ca;?~"^e^^ , as a; — > oo; (a;) Ca;~?, as a; —> —oo 

The asymptotics of fpix) = <y?+(— a;) is obvious. 

The boundary behavior of the process and Y is given by the following lemma: 

Lemma 4.13. Let h{x) = ci(p+(x) + C2(p~{x). Then both boundaries are natural for the process X^ 
(the same for the process Yt). 

Proof. Let's prove the result for the right boundary — oo. The result for Z)^ = — oo follows by 
symmetry. 

Let ci > 0. From lemma 3.8 and equations (B.2) we find that the speed measure and the scale function 
of X^ have the following asymptotics as a; ^ oo: 

mxh(x) = h^{x)mx{x) ~ Cx^'e^^', s'xh{x) = h-^{x)s'x{x) ~ Cx-^^e-^"^', 

where q = f — 1. Now using lemma 2.7 one can check that = oo is a natural boundary. The case 
ci = and C2 > can be analyzed similarly. □ 

Thus we see that the processes Yj associated with the Omstein-Uhlenbeck process behave similar 
to processes associated with Brownian motion (quadratic volatility family): these are conservative pro- 
cesses. The next two examples of processes illustrate different boundary behavior: the family associated 
with the CIR process has one killing (exit) boundary, while the family associated with the Jacobi process 
has both killing (exit) boundaries. 

4.4. CIR processes - confluent hypergeometric family of driftless processes. Let Xt be the CIR pro- 
cess 

dXt = (a - hXt)dt + G^/XtdWt, 
considered in section B.2. We will use the notations from section B.2: 

a = ^ — 1 and = 

Functions <p+ and ip~ for the CIR process are solution to the ODE 

(4.18) .a^^^^-{a-hx)^^p^(x). 

Making affine change of variables y = Ox and dividing this equation by b, we reduce equation (4.18) to 
the Kummer differential equation (A.ll), thus and ip'^ satisfy 

(4.19) ^+{x)^M{^,a + l,ex), 

(4.20) ^-{x) = U{^,a+l,ex). 

Using formulas (A.13) and (A.15) we find the asymptotics of (Pp{x) and (Pp{x): 

(4.21) (/?+(x) ~ 1, asx^O; <^+(x) ~ Ce**="x^-"-\ as x ^ oo; 

(4.22) (f~ (x) ^ Cx"" , as X ^ 0; ip~ (x) ^ Cx~^ , as x ^ oo. 

The following lemma describes the boundary behavior of the process X'^ (and thus of the process 

Yt = y(x")). 

Lemma 4.14. Let h{x) = Ci(p+(x) + C2(p~{x), where both Ci are positive. Then = oo is a natural 
boundary of the process X^, while = Q is a killing boundary if a e (0, 1) and it is an exit boundary 
if a > I. 
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Proof. Using formulas (4.21),(B.9) and lemma 3.8, we find that the asymptotics of the speed measure 
and scale function of X'* are 



mxh{x) = h {x)mx{x) 

and 



Cx^h " ^gfla;^ as X ^ OO 

Cx~", as a; — > 0, 



s' (x)-h-^(x)s' asar^oo 

One can check using Feller's theorem 2.7, that for a. E (0, 1) we have a regular boundary, but since the 
/i-transform introduces nonzero killing measure at the boundaries given by equation (3.12), we have a 
killing boundary. If a > 1 we have an exit boundary. □ 

Note that since the left boundary Z)^ = is either a kilHng or an exit boundary, the process X'* is not 
a conservative process. The same is true for = Y{Xf). 

Theorem 4.15. The family ofdriftless processes related to a CIR process by a stochastic transforma- 
tion is characterized by their volatility functions as follows: 



(4.23) (^y{Y{x)) = 



{ciip'j,{x) + C2(pp (a;))2' 



ci(/9^(a;) + C2ipp{x) 

Definition 4.16. We will call this family of driftless processes the confluent hypergeometric family and 
denote it by 

m{ClR) = {Yf : y - CIR process}. 
4.5. Jacobi processes - liypergeometric family of driftless processes. Let Xt be the Jacobi process 

dXt = {a - bXt)dt + a^X,{A-X,)dWt^ 
described in section B.3. Recall the notations from section B.3. We defined parameters 

26 2a 2a 
a = ^ ^ - 1 and /3 = -^r— - 1. 

We assume that a > and (3 > 0, thus both boundaries are inaccessible (see section B.3). 
9+ and ip- 



Functions and (p for the Jacobi process are solutions to the ODE 



By the affine change of variables x = Ay this equation is reduced to the hypergeometric differential 
equation (A.2), thus using equations (A.8) we find that functions and (p~ for the Jacobi process are 
given by 

(4.24) ip+{x) = 2Fi{ai,a2;Pi;x/A) 

(4.25) ip-{x) = 2-Fi(Q!i, as; ai + as + 1 - 1 - x/A) 
where the parameters satisfy the following system of equations: 

'ai + a2 + l = fl 

(4.26) I aia2 = f§ 

The asymptotics of iPp{x) and iPp{x) are 

(4.27) ~ 1, asa; 0, ~ C(A - a;)"", as a; ^ A, 

(4.28) ^p{^) ^ Cx~^, as a; ^ 0, <pp{x) ^ 1, as a; — > A. 
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Let h{x) = ci</?+(a;) + C2fp (x), where both Cj are positive. The boundary behavior of the h- 
transformed process Xj^ is similar to the case of CIR process at x = 0: 

Lemma 4.17. = is a killing boundary for the process if (3 G (0, 1) md it is an exit boundary 
if 13 > 1. The same is true for D"^ = A by changing (3 a. 

We see that in the case of the Jacobi process both boundaries are either killing or exit boundaries for 
Xf, thus Xf and are not conservative processes. 

Theorem 4.18. The family ofdriftless processes Yt related to a Jacobi process by a stochastic transfor- 
mation is characterized by their volatility functions as follows: 



(4.29) <yY{Y{x)) = C^/x{A - x) 



X 



-a-l 



{A-x)-f^-^ 



{cnp'^{x) + C2ipp 

^, ^ C3ip+{x) + C4,ip-ix) 

Y{x) — 



CHpj{x) + C2iPp {x)' 

Definition 4.19. We call this family of driftless processes the the hypergeometric family and denote it by 

SDT(Jacobi) = {Yt : Y ^ Jacobi process}. 

5. Classification of driftless diffusion processes 

In the introduction to section 3 we discuss the transformations of diffusion processes that preserve 
the solvabiUty property. As we saw these transformations consist of change of variables and Doob's 
h-transform, which also can be considered as a change of variables and gauge transformation of the 
Markov generator. In this section we will focus on Markov generators. Let's ask a question: how can one 
transform operator £ and solutions to the eigenfunction equation? 

Let C be the second order differential operator 



(5.1) C = a{y)^ + b{y)- 



If a{y) is positive on some interval D operator C can be considered as a generator of a diffusion process, 
thus we will call a{y) the volatility coefficient and b{y) the drift coefficient of operator £. 
Let's consider the following three types of transformations of the operator C: 

(i) Change of variables x = x{y): the solution of the eigenfunction equation Cf = pf is mapped 
into f{y) f{y{x)) and 

(5.2) £ ^ Ty^^C = a{y){x'{y)f^ + {a{y)x"{y) + b{y)x'{y))^ = 

(5.3) aiy)ix'iy)r^ + {Cx){y)^, 

where y = y{x). 

(ii) Gauge transformation: f{y) t-^ f{y)/h{y) and 

d"^ { h'(v]\ d 1 

£^n£= ICh = a{y) ^ + (h{y) + 2a{y) -^j - + - {a{y)h"{y) + b{y)h'{y)) = 

(5.4) = aiy)^ + (%) + 2a(,)^) |- + J-^(£/.)(,). 

Notice that gauge transformation actually consists of two transformations: right multiplication 
of £ by /i and left multipUcation by l//i. 

(iii) Left multipUcation by a function 7^(0;) (which does not affect /(a;)): 

(5.5) £ - T^2£ = ^^(y)a{y)^+'y^{y)b{y)^. 
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In the case when a{y) = ^^^(j/) and £ is a generator of a Markov diffusion (Yj.P), we can give a 
probabilistic interpretation to the first and second types of transformations described above: Ty^^ is just 
the usual change of variables formula for the stochastic process Yt, which describes the dynamics of the 
process Xt = X{Yt), thus it is just an analog of the Ito formula written in the language of ODEs. The 
Gauge transformation has a probabilistic meaning if h{Yt) can be considered as a measure change density 
(thus h{Yt) is a local martingale and Ch = 0), or when /i is a p-excessive fiinction (Ch = ph) - then 
the gauge transformation \Ch — p is the Doob's /i-transform discussed in section 3.2. The transformed 
generator C in this case describes the dynamics of the process Yj under the new measure Q, defined by 
dQt = h{Yt)<IPt- Note that both of these transformations preserve the form of backward Kolmogorov 
equation: 

|/(t,y) = £/(i,t/). 

The last transformation T^2 does not preserve the form of backward Kolmogorov equation, thus in 
general it has no immediate probabilistic meaning, except when 7^(0;) = c is constant - then T^2 is 
equivalent to scaling of time ^t' . 

As we have seen in section 3, the OU,CIR and Jacobi families of processes are solvable because they 
can be reduced to some simple solvable process. In other words, for these processes, the eigenfunction 
equation 

(5.6) CYf{y) = laUy)^0^=Pm 

can be reduced by a gauge transformation (change of measure) and a change of variables to a hyper- 
geometric or a confluent hypergeometric equation, thus giving us eigenfunctions ipn{x), generalized 
eigenfunctions ip~^,(fi^, and a possibility to compute the transitional probability density as 

(5.7) pY{t,yQ,yi) = -r- — px{t,xo,xi) = — — -r- — r V e~^''"'^'''Vn(a;o)V'n(a;i). 

h{xo)h{xi) h[xo)h[xi) ^ 

It is known that OU, CIR and Jacobi processes are the only diffusions associated with a system of or- 
thogonal polynomials (see [11]), thus corresponding families of driftless processes are the only ones that 
can have a probability density of the form (5.7) where the orthogonal basis {ipn}n>o is given by orthog- 
onal polynomials. However we might hope to find new famiUes of solvable processes if we generalize 
the definition of solvability. 

Note that the fact that we can compute solutions of equation (5.6) gives us a ready expression for the 
Green function through the formula (2.19): 



(5.8) GYiX,yo,yi) = 



wx^vtiyi)vxiyo), yi < yo- 



Since the Green function is a Laplace transform of (t, yo, yi) one could hope to find the probability 
kernel through the inverse Laplace transform of G'i'(A, yo, Thus in this section we will use the 
following definition of solvability: 

Deiuiition 5.1. The one dimensional diffusion process Yt on the interval Dy is called solvable, if its 
Green function can be computed in terms of (scaled confluent) hypergeometric functions. 

In other words, the process Yt is solvable if there exist a A-independent change of variables y — y{z) 
and a (possibly A-dependent) fimction h{z, A), such that all solutions to the equation 

(5.9) jCyfiy) = Xf{y) 

are of the form h{z{y), X)F{z{y)), where F is either a hypergeometric function 2-^1(0, b; c; z), or a 
scaled confluent hypergeometric function M(a, b, wz) (with parameters depending on A). 

Remark 5.2. Later we will see that the requirement that the change of variables is independent of A is 
necessary, because otherwise every diffusion process is solvable (see remark (5.14)). 
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5.1. Liouville transformations and Bose invariants. Consider the linear second order differential op- 
erator 

(5.10) Cy = a{y)^ + biy)-^. 
Making the gauge transformation with gauge factor h 

we can remove the "drift" coefficient and thereby arrive at the operator 

1 

£y ^ -Cyh = a(y)^ + a{y)I{y). 
Multiplying this operator by a{y)~^ we arrive at symmetric operator 

(5.11) ^Lyh ^ a-'^Cyh = jC<^=^+ I{y), 
where the potential term is given by: 

26(y)a(y)' - 2a(y)6(y)' - biy)' 



J, . (h'{y)\ (h'{y)\ 



Definition 5.3. C is called the canonical form of the operator Cy . 

The form of operator is clearly invariant with respect to any of the three types of transformations 
described above. Moreover, the following lennma is true: 

Lemma 5.4. The canonical form of the operator given by (5.11) is invariant under any two transfor- 
mations of {Ty^z , Th , }. 

Proof. This lemma is proved by checking that every combination of two transformations cannot change 
the canonical form of the operator (5.11). For example, suppose we are free to use Ty^^ and Tf^, but 
not T^2. Using formula (5.4) we find that by applying Th we have nonzero drift given by 2h'/h. We 
can't remove this drift by some change of variables Ty^x, since by formula (5.2) it will add a nontrivial 
volatihty term x'{y)^, which in turn can not be removed by any gauge transformation T^. As another 
example let's assume that we can use Ty^^ and T^2, but not Th. By formula (5.5) we see that T^2 adds 
nontrivial volatility, which can be removed by Ty^x, but formula (5.2) tells us that this in tum will add 
nonzero drift ■^^{x){£x){y), which can't be removed by any T^2. The last combination T^2 and Th can 
be checked in exactly the same way. □ 

Remark 5.5. Note that to bring Cy to canonical form we used two types of transformations: a gauge 
transformation and a change of variables. But by using other choices of two transformations we could 
bring Cy to a different canonical form. Thus when we talk about canonical form we need to specify 
with respect to which two types of transformations this form is invariant (in this section we use only two 
types of canonical forms: one described above, and the second obtained by a gauge transformation and a 
change of variables only). 

Definition 5.6. Function I{y) is called the Bose invariant of operator (5.10) (with respect to transforma- 
tions Th and T^2). 

As we proved, function I{y) is invariant with respect to any two transformations of {Ty^x, Th, T^^ }. 
However it is possible to change the potential I{y) by applying all three types of transformations. The 
idea is to apply first a change of variables transformation, then remove the "drift" by a gauge transfor- 
mation, after which we divide by the volatiUty coefficient to obtain a new canonical form. The details 
are: 

(i) A change of the independent variable y = y{x) (see equation (5.2)) changes the operator (5.11) 
into 
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(ii) Multiplying the operator by 7 (a;) — {y'{x)) we arrive at 

(iii) Applying the gauge transformation with gauge factor h = \/y'{x) (see equation (5.4)) brings us 
to the operator in the following canonical form 

~ 'q^ + J{x), 

where the potential term is transformed as 
(5.13) J(x) = \{y, x) + {i)\x)fl{i,{x)\ 

and {y, cc} is the Schwarzian derivative ofy with respect to x: 



y'{x)J 2\y'{x) 

As we have seen the above transformation changes the canonical form of the operator by applying all 
three types of transformations. It is called a Liouville transformation. 

Note that the order of the different steps in Liouville transformation does not matter, since all the three 
transformations {Ty^^, Th, Tlys} commute. The other important idea is that we are free to choose the 
first transformation, but the other two are uniquely defined by the first one (see lemma (5.4)). We will 
use this fact in the proof of the main theorem in the next section. 

As the first application of the canonical forms and Bose invariants we will prove the following leirnna, 
which gives a convenient criterion to check whether two processes are related by a stochastic transfor- 
mation: 

Lemma 5.7. Let be a diffusion process with Markov generator 

^ ^ ,/ ^df(x) 1 .d'^fix) 

With this diffusion we will associate a function Jx{z) = Ix{x{z)), where 

(5.14) Ix{x) = \ (aix)a"ix) - '-(a^x))^ + 2 [^^^ _ - ^]) . 

and the change of variables x{z) is defined through x' {z) = a{x). Then the function Ix{z) is an 
invariant of the diffusion Xf in the following sense: 

(i) Yt = {Y{Xt),P) ifandonlyifJriz) = J x{z) for all z. 

(ii) X^ ~ [Xt , P'') is an p-excessive transform of Xt if and only if Jx>^ {z) = Jx [z) — p. 

(iii) X ^ Y if and only if Jy{z) = Jx{z) — p, where p is the same as in stochastic transformation 
{p, h, Y} which relates X and Y. 

Proof. Applying the gauge transformation with gauge factor hi (x) = ^/W(x) changes Cx as follows: 

(5.15) Cx ^ ^rx/ii = \^xix)^ + h{x), 
where by formula (5.12) the potential term is equal 

ibx^xy^jx) _ ^ _ h\{x) 

^x{x) ^ a\{x) 

Changing variables x\z^ = ^'^xix) we arrive at 

(9^ r"{z\ d 



(5.16) /i(x) = i 
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and at last making the gauge transformation with gauge factor /12 (-2) = \/x'{z) we arrive at the canonical 

form 

(5.18) £^ = ^£,/i2 = ^ + Jx{z) 
where the Bose invariant (potential term) is 

(5.19) Jx{z) = ]^{x,z} + h{x). 

By direct computations we check that {x,z} = \(Jx{x)ax{x) — \(T\{x), which gives us formula (5.14). 

Now we need to prove that Jx{z) is invariant under stochastic transformations up to an additive 
constant. A change of variables does not change the form of Jx{z). A change of measure (p-excessive 
transform) changes the generator of X 

Cx ^ \^x^ ~ P- 

Since Jx{z) is invariant under gauge transformations Cx 1— > j^^xh we see that Jy{z) differs from 
Jx {z) by a constant —p, which ends the proof. □ 

Example 5.8. As we showed in the previous section the quadratic volatility family with volatiUty func- 
tions of the form 

(5.20) cry(y) = a2y^ + ai2/ + ao 

is related to a Brownian motion Xt = Wt by a stochastic transformation. It is very easy to prove this fact 
using the above lemma: 

J\Y = 0, Jy = const, 

thus Ytr^Wt. 

Example 5.9. As another apphcation of lemma 5.7 let's show that a Bessel process 

dXt = adt + ^/XtdWt 
is related to CEV processes (constant-elasticity-of-variance) 

dYt = YfdWt. 

Note that in subsection 4.2 we constructed these stochastic transformations explicitly. 
Function Ix {x) is given by 



a 


a"- 




X 


X 


) 



^(_2a2 + 2a-3/8) 



After making the change of variables dx = ax{x)dz = \fxdz (thus x = -2^/4) we arrive at 

Jx{z) = z-2(-2a^ + 2a - 3/8). 

Similarly for the process Yj 

dy = (JY{y)dz = y'^dz ^ y = {{1 - 9)z)^ 

and 



Mz) = -Ae{e-i)y^'-^-'-6^y'^--\ = '-z-^ 2 



4 {i-ey 

One can check that when 9=1— 2(2l-i) ^^^^ Jy{.z) = Jx{z), thus processes Yt should be 
related to the Bessel process by a stochastic transformation {p, h, Y} with p = (see section 4.2 where 
we construct this transformation explicitly). 
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The following diagram illustrates the importance of a Liouville transformation (in particular the T,^2 
transformation). Assume that we start with a diffusion process Xi . By stochastic transformations we 
can construct the family S!Jt(Xi) of driftless processes. By a gauge transformation and by a change of 
variables we map the generator of Xi into the corresponding canonical form (with the Bose invariant 
given by (5.14)). By applying a Liouville transformation we change the potential (Bose invariant) of the 
canonical form and then by a change of variables and a gauge transformation we can map it back into 
the generator of some diffusion process X2. Then as before we construct a family of driftless processes 
d}l{X2). Note that Tl{X2) and Tl{Xi) are not related by a stochastic canonical transformation, since 
otherwise they would have the same Bose invariants. Then the process can be continued. Thus we have 
a family of Bose invariants related by Liouville transformations, each of these invariants (potentials in 
canonical form) gives rise to a new family of driftless processes, which are not related by a stochastic 
transformation to the previous families (otherwise their Bose invariants would coincide). 

The following diagram summarizes the above ideas: 



I 

stochastic 
transformations 

; 

Xi 

t 



7^ 



Liouville 
transformation 



9K(X2) 

I 

stochastic 
transformations 

I 

t 



7^ 



Liouville 
transformation 



m{Xs) 

I 

stochastic 
transformations 

^3 
t 

Jx^z) 



5.2. Classification: Main theorems. 

Theorem 5.10. First classification theorem. A driftless process is solvable in the sense of definition 
5.1 if and only if its volatility function is of the following form: 



(5.21, = ',yiyi^)) = c^ \,^,JtU^^)r V WY 

where the change of variables is given by 

(5.22) y = y(a;) = — — — — — , CiC4-C3C2^0. 

ciFi (a;) +02^2(0;) 

In the case A{x) = x: 

(i) R{x) € P2, such that R{x) ^ in (0, 00) 

(ii) Fi and F2 are functions M{a, b, wx) and U (a, b, wx) 

(iii) W{x) is a Wronskian of the scaled Kummer differential equation (equal to s'{x) in (B.9)). 
and in the case A{x) = x{l — x) 

(i) R{x) g P2, such that R{x) ^ in (0, 1) 

(ii) Fi and F2 are two linearly independent solutions to the hyper geometric equation given by (A.3). 

(iii) W{x) is a Wronskian of the hypergeometric differential equation (equal to s'{x) in (B.17)). 

Before we prove this theorem we need to establish some auxiliary results. 
Lemma 5.11. The Bose invariant for the hypergeometric equation 

x{l - x) + (-y - (1 + a + (3)x) ^ - am^) = 0- 

is given by 
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where 

(5.24) Q{x) = (1 - (a - /S)^),^^ + (27(0 + /3 - 1) - 4a^)a; + 7(2 - 7). 
77ie Bose invariant for the scaled confluent hypergeometric equation 

(5.25) honflix) = 
w/zere 

(5.26) Q{x) = -w^x^ + 2w{a - 2h)x + a{2 - a). 

In both cases by varying parameters we can obtain any second order polynomial Q. 
Proof. For the second order differential equation 

(5.27) aix)^^ + Kx)^ + c{x)f{x) = 

the Bose invariant is given by formula (5.12) and is equal to 

_ 2h{x)a{xy - 2a{x)b{xy - b{xf c{x) 

_ 2b{x)a{x)' - 2a{x)b{x)' - b{xf + 4a(x)c(x) 
4a(a;)2 



□ 



Corollary 5.12. Let T{x) G P2 be an arbitrary second order polynomial and A{x) G {x,x{l — x)}. 
The solutions to equation 

can be obtained in the form F{x) / yJW{x), where W{x) is the Wronskian given by exp(— ^^dy) 
and function F{x) is a solution to the hypergeometric equation in the case A{x) = x{l — x) or to the 
scaled confluent hypergeometric equation in the case A{x) = x. 



(5.28) ^{y,x} = J{x) 



Theorem 5.13. (Schwarz) The general solution to the equation 
(5.28) 
has the form 

(5.29) y{x) = fM, 

Fi{x) 

where Fi and F2 are arbitrary linearly independent solutions of equation 

(5.30) + J{x)F{x) = 0. 
Proof. Let's introduce the new variable F(x) = J^— . Then we have 

y'(-) = lF^ -iog,'(x) = -2iog(^(x)) -^ = -25^. 

Thus the Schwarzian derivative {y, z} is equal to: 

,'y"{x)y 1 fy"{x)Y ^F"{x) 
{y,z} ' 



y'{x)J 2\y'{x)J F{x) 
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and equation (5.28) is ti-ansformed into equation 

F"{x) + J{x)F{x) = 

for function F{x), thus 

1 



y'ix) 



where Fi{x) is an arbitrary solution of F"{x) + J{x)F{x) = 0. 

Note that if Fi and F2 are solutions to second order Unear ODE aF" + bF' + cF = 0, then 

d fF2{x)\ F^{x)F,{x) - F2{x)Fi{x) Wf,,fA^.) 



dx \Fi{x) J F^{x) Ff{x) ' 

and since the Wronskian of equation (5.30) is constant we can obtain the above expression for y{x). □ 

Remark 5.14. Note that theorem 5.13 and equation (5.13) tells us that for any potential there exists a 
change of variables z{x), which maps equation F" = into equation F"{x) + J{x)F{x) = 0, thus any 
two canonical forms can be related by some Liouville transformation. That is why in the definition (5.1) 
we require the change of variables y{x) to be independent of A. 

Proof of the first classification tlieorem: 



Proof. By definition 5.1 the process is solvable if for aU A we can reduce equation 



to the (scaled confluent) hypergeometric equation by a Liouville transformation, with a change of vari- 
ables y = y{x) independent of A. The Liouville transformation consists of three parts: a change of 
variables, a multipUcation by function and a gauge transformation and all these transformations com- 
mute. For example we could first apply a multiplication by function transformation, and then use only 
change of variables and gauge transformations. Thus we can reformulate the problem as follows: find all 
functions (Jviy), such that there exist a function 7(2/), such that the equation 



(5.31) 2^\y)CYf{y) = ^Y{y)i\y)^^ = '^^i\y)f{y) 



dy'^ 

can be reduced to the (confluent) hypergeometric equation by a change of variables (independent of A) 
and a gauge transformation. The Bose invariant of equation (5.31) is given by 

(5.32) 7(ar) = i{y,x}-2A7'(y(x)), 
with the change of variables given by 

(5.33) £ = ar(y)7(y)- 

Note that 7(2/) must be independent of A since cry (y) and y'{x) are independent of A. 

By lemma 5.4 we know that equation (5.3 1) can be reduced to the (confluent) hypergeometric equation 
by a change of variables and a gauge transformation if and only if the corresponding Bose invariants are 
equal, that is 

I{x) = Ihyp{x), or I{x) = Iconfl{x), 

thus, combining equations (5.32) with (5.23) and (5.25), we have 

]{y,x}-2X^\y{x))- ^^""^ 



where A{x) G {x, .t(1 — x)} and Q{x) = Q{x, A) is some second order polynomial in x with parameters 
depending on A. Note that {y, x} and j'^{y{x)) are independent of A, thus there exist two polynomials 
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T{x),R{x) e P2, independent of A, such that 

{Q{x) =T{x) -2XR{x), 
7'(y(^)) = i^, i?(x)>0. 
The last equation in the system (5.34) gives us the function 7(2/(2;)): 



(5.35) jivix)) 



A{x) 



By theorem 5.13, the solutions to equation ^{y, x} = are given by 



(5.36) y{x) 



cshix) + 04/2(2;) 



cih{x) + 02/2(2;)' 
where /i and /2 are linearly independent solutions to equation 

(5.37) fix) + = 0. 

Now we can use corollary (5. 12), which teUs us that aU the solutions to equation (5.37) can be found in 
the form F{x)/^/W(x), thus 

(5.38) vix) = ^3^i(^) + ^^g(^) 

' CiFi{x) + C2F2{x) 

where F-i and F^ are confluent hypergeometric (A{x) = x) or hypergeometric iA{x) = x{l — x)) 
functions. 

Now we are ready to find volatihty function crY(y). From the equation (5.33) we find that 

(5.39) aY{Y{x)) = y'{x)- ^ 



7(y(a;))' 

The derivative y'{x) can be computed as 

CW{x) 



y'{x) = 



{ciFi{x) + C2F2{x))^' 

thus 

ar{Y[x)) - y (-)^(^ - ^ i^c,F,{x) + C2F2{x)Y 

which completes the proof. □ 

Definition 5.15. We call the family of driftless processes with volatility function given by equation (5.21) 
a hypergeometric R-family in the case A{x) = x{l — x) and a confluent hypergeometric R-family in the 
case A{x) = x. 

We see that in the case R{x) = A{x) we recover the hypergeometric and confluent hypergeometric 
famihes, which correspond to 9Jl(Jacobi) and 9Jl(CIR). In the case R{x) ^ A{x) we obtain new families 
of processes. The next theorem shows that the (confluent) hypergeometric R-family can be obtained by 
stochastic transformations described from some diffusion process (in the same way as Jacobi and CIR 
families are generated by a single diffusion process). 

Theorem 5.16. Second classification theorem. Let R{x) E P2 be a second degree polynomial in x. 

(i) The confluent hypergeometric case: Assume that R{x) has no zeros in (0, 00). Let Xt = Xl^ be 
the diffusion process with dynamics 

(5^40, + 
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Then the confluent hypergeometric R-family coincides with dyt{X^) and thus can be obtained 
from by stochastic transformations. In the particular case R{x) = A{x) = x we have 
is a CIR process and we obtain the CIR family defined in (4.4). 
(ii) The hypergeometric case: Assume that R{x) has no zeros in (0, 1). Let Xf = X^ be the 
diffusion process with dynamics: 

(5.41, .X. = (»..X.)^^* + ^^«. 

Then the hypergeometric R-family coincides with DJl{X^) and thus can be obtained from X^ 
by stochastic transformations. In the particular case R{x) = A{x) = x{l — x) we have Xp is 
a Jacobi process and we obtain the Jacobi family defined in (4.5). 

Proof. The generator of Xt is given by 

(5.42) ^- = (« + ^-)4)a. + 2Wa^- 

One way to prove this theorem is to check that the corresponding Bose invariants coincide. However we 
prove this theorem by applying stochastic transformations to the process Xt and showing that we can 
cover all of the processes with volatiUty functions given by equation (5.21). 

First we need to find two linearly independent solutions to the "eigenfunction" equation 

CxV = pip- 

This equation is equivalent to 

, .difiix) .d'^uDix) „ R(x) , . 

2(a + 6.)^+A(.)^ = 2p^^(.). 

By dividing both sides by A(x) and making gauge transformation with gauge factor h(x) = \fW(x), 
F = (p/h,we arrive at the equation in canonical form: 

d^F(x) Q{x)^2aR{x) 
dx^ ^ A^x) ^ ' 

By corollary 5.12 this equation is solved in terms of hypergeometric functions, thus (pi{x) = g{x)Fi(x), 
where Fj are (confluent) hypergeometric functions. Thus 

Yf^ _ Ci(/?i(x) + C2<f2{x) _ CiFi(x) + C2F2{x) 
C3^l{x) + C4,<P2{x) C3Fi(x) + C4,F2{x)' 

and ayiu) is computed as 



ariy) = Y'{x)4^ = C^/Aixj- 



W{x) A{x) 



^/RXx) " ^ '{crF^(x)+C2F2{x)Y\l R{xy 
which ends the proof. □ 
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Appendix A. Hypergeometric Functions 

In this section, we review basic notions about hypergeometric functions. A good collection of facts 
and formulas can be found in [8] and [1]. 

A.l. Hypergeometric function. The hypergeometric function 2Fi{a,(5;^;z) is defined through its 
Taylor expansion: 



(A.1) 2Fi(a,/3;7;z) = ^ 



(a)„(/?)„ z" 



2-Fi(q;, /3; 7; z) is a solution to the hypergeometric differential equation 

(A.2) z{l - z)F"{z) + (7 - (1 + a + P)z)F'{z) - a0F{z) = 0. 

This differential equation has three (regular) singular points: 0, 1, 00. The exponents at 2; = are 0,1 — 7 
and at z = 1 are 0, 7 — a — /?. 

Two Unearly independent solution in the neighborhood of z = are given by: 

(A.3) wi=2Fi{a,p;r,z), 

(A.4) W2 = z^-^2Fi{a - 7 + 1, /3 - 7 + 1; 2 - 7; z), 

and in the neighborhood of z = 1 

(A.5) wi = 2Fi{a,l3;a + f3 + l-rA-z), 

(A.6) «;2 = (l-^)^-"-^2-Fi (7-/3,7- a,7- a- /3+ 1,1-^)- 

The derivative of the hypergeometric function is: 

(A.7) 2Fiia,l3;r,z) = ^2Fiia + l,p+ln+l,z). 

7 

Increasing and decreasing solutions of the hypergeometric equation, which in the case a > 0, /? > 
0, 7 > and 7<Q! + /3+lare given by: 

(A.8) ip+ix) = 2Fi{a,P;r,z), 

(A.9) ip-{x)=2Fi{a,0;a + P+l-rA-z). 

A.2. Confluent hypergeometric function. The confluent hypergeometric function iFi{a,b,z) (also 
denoted by M{a, b, z) or $(a, 6, z)) can be defined through its Taylors expansion 



(a)n 
{b)n n\ ■ 



(A.10) M{a,b,z) = Y, 

n=0 

Function M(a, b, z) is a solution to the Kummer differential equation 

(A.ll) zF"{z) + {b- z)F'{z) - aF{z) = 

It has two singular points: 0, 00. is a regular singular point with the exponents 0,1 — 6. 
Two Unearly independent solutions to the Kummer differential equation are given by: 

(A.12) wi = M{a;b;z), W2 = z^-^M{l + a - b;2 - b; z). 

The increasing and decreasing solutions are given by: 

(A.13) ip+{x) = M{a\b\z), 

/A -/ N rr/ I, N ^ / M{a,b,z) i_f, M(l + a - 6, 2 - 6, z) " 

(A.14) ip {x) = U{a,b,z)= ^-—[——^ T^^-^ 



sin(7r6) Vr(l + a - b)T{b) T{a)T{2 - b) 
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The asymptotics of M and f7 as |2;| — > oo > 0) is: 

m 



(A.15) ^+{z) = M{a, b, z) = ^e^z"-\l + 0{\z\-')), 

(A.16) ^-{z) = U{a, b, z) = z--{l + 0{\z\-')), 

and the derivative of the confluent hypergeometric function can be computed as 

(A.17) M'{a,b,z) = '^M{a+l,b+l,z). 

Appendix B. Ornstein-Uhlenbeck.CIR, and Jacobi processes 

In this section we present some facts about the following three diffurion processes: OU, CIR and Ja- 
cobi. These processes enjoy a lot of interesting properties: first of all , these processes are associated with 
a family of orthogonal polynomials, which means that these polynomials form a complete set of eigen- 
functions of the generator C Furthermore, it can be proved that these are the only diffusion processes 
which have this property (see [11]). This property of generator £ allows us to express the probability 
semigroup P{t) as an orthogonal expansion in these polynomials, thus all of these processes are solvable 
(moreover, there are explicit formulas for OU and CIR). 

In this section we just briefly present the formulas for the generator, speed measure and scale func- 
tion, describe boundary behavior, present eigenfunctions and eigenvalues of the generators and series 
expansion for the probability density (along with an expUcit formula if possible). 

B.l. Omstein-Uhlenbeck process. 

• Generator: 

where b > 0. 

• Domain: D = (— oc'.+oc). 

• Speed measure and scale function: 

(B.2) m{x)= —L=cxp(^-^{x-^y^ , s'{x) = exp(^^{x-^y^ 

• Boundary conditions: Both = — oo and D"^ = +oo are natural boundaries for all choices of 
parameters. 

• Probability fimction: 

iB.^^^^\t,xo,xi)m{xi) = , .exp ^- ^ . 

• Spectrum of the generator: 

(B.4) A„ = -bn. 

• Eigenfunctions of the generator: 

(B.5) M^) 

where Hn{x) are Hermite polynomials. The three term recurrence relation is: 
(B.6) - 2xHn{x) + 2nHr,-i{x) = 0. 

• Orthogonality relation 

ijjn{x)tpmix)m{x)dx = 2"n!(5„ 
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• Eigenfunction expansion of the probability function: 
(B.7) p^°^\t,xo,xi) = J2 '"^HMHniyi), 



2"n! 

ri=0 



where = ^ ^{xi - f). 

B.2. CIR process. 

• Generator 

(B.8) ''=^'-'-^i + r'-&^- 

• Domain D = [0, +00) 

• Speed measure and scale function: 

(B.9) m{x) = J"^\^ x"e-'-, s'{x) = x—'e^'^ , 

where Q!=^-land6'=f|. 

• Boundary conditions: D'^ = +00 is a natural boundary for all choices of parameters and 

exit, if a < — 1 
regular, if — 1 < a < 
entrance, if < a 



(B.IO) = < 

• Probability function: 



(bt\ 2"-^ / 
j exp [-Ct(a;oe-" + xi)] (204^X0X16" 

wherect = — 26/(c7^(e~''* — 1)) and/a is the modified Bessel fimction of the first kind (see 

• Spectrum of the generator: 

(B.12) A„ = -bn. 

• Eigenfunctions of the generator: 

(B.13) M^') = KiOx), 

where (y) are Laguerre polynomials of order a with the three term recurrence relation: 
(B.14) (n + l)LZ+,{y) - (2n + a + 1 - + (n + a)L^_,{y) = 0. 

• Orthogonality relation: 

1pn{x)lprn{x)m{x)dx = '^7^(5„m- 

D nl 

• Eigenfunction expansion of the probability function: 

(B.15) p^^'^\t,xo,x,) = Y^e-'"^'-^^LZ{exo)LZ{ex^). 

„=o ^'^ + 

B.3. Jacob! process. 

• Generator 

(B.16) C= (a-bx)-^ + -a'^x(A-x)-^. 

ax 2 ctx^ 

• Domain L» = [0,A] 

• Speed measure and scale function: 

= ^.7;t7^I?I^:^4";^. Ax) = x-f^-\A- x)—\ 



where a=||-^-land/3=^-l. 
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(B.18) 



• Boundary behavior for the Jacobi process is the same as for CIR process at the left boimdary: 

'exit, if /? < -1 
regular, if - 1 < /3 < 
^ entrance, if < /3 

The same classification applies to right boundary, we only need to replace f3 by a. Notice that 
in the case when a > 0, 6 > and | < A (which means that mean-reverting level lies in the 
interval (0, A)), we have a > —1 and /3 > — 1 and thus both boundaries are not exit. 

• Spectrum of the generator: 

(B.19) A„ = -^n(n-l+^). 

• Eigenfunctions of the generator: 

(B.20) = PjC'^'Hy). 

where y = ( ^ — 1) and Pi°''^^ (y) are Jacobi polynomials with the three term recurrence relation: 

yp("'0)(y) - 2(n + l)(n + Q + /3+l) M)r^^ , 



(2n + a + /3)(2n + a + /3 + 2) " 
, 2(n + a)(n + /3) («,/?) 

^{2n + a + P){2n+a + P+l) ""^ 



Orthogonality relation 



(a + p + 2)„_i(2n + a + p + Ijn! 



ID 

• Eigenfunction expansion of the probability function: 

oo —\t 

(B.21) p(J-*)(i,xo,a;i) = ^ ^-^Pi^'''\yo)P^"^^\y^), 



n=0 



where yi = - 1) 
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